Relaxation times in an open interacting two-qubit system 



OO 
O 

o 

(N 
Oh' 

OO 



Y. Dubi and M. Di Ventra 
Department of Physics, University of California San Diego, La Jolla, California 92093-0319, USA 

In a two-qubit system the coupling with an environment affects considerably the entanglement 
dynamics, and usually leads to the loss of entanglement within a finite time. Since entanglement 
is a key feature in the application of such systems to quantum information processing, it is highly 
desirable to find a way to prolonging its lifetime. We present a simple model of an interacting 
two-qubit system in the presence of a thermal Markovian environment. The qubits are modeled as 
interacting spin-i particles in a magnetic field and the environment is limited to inducing single 
spin-fiip events. A simple scheme allows us to calculate the relaxation rates for all processes. We 
show that the relaxation dynamics of the most entangled state exhibit critical slowing down as a 
function of the magnetic field, where the relaxation rate changes from exponentially small values to 
finite values in the zero-temperature limit. We study the effect of temperature and magnetic field on 
all the other relaxation rates and find that they exhibit unusual properties, such as non-monotonic 
dependence on temperature and a discontinuity as a function of magnetic field. In addition, a 
simple scheme to include non-Markovian effects is presented and applied to the two-qubit model. 
We find that the relaxation rates exhibit a sharp, cusp-like resonant structure as a function of the 
environment memory-time, and that for long memory-times all the different relaxation rates merge 
into a single one. 

PACS numbers: 03.67.Mn, 03.65.Yz, 03.65.Ud, 42.50.Lc 
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I. INTRODUCTION 

When considering application of quantum information 
processing, two main ingredients must be considered - 
entanglement and decoherence. While entanglement, or 
non-local coherence, plays a key role in qubit operations 
decoherence sets the limit to which such operations 
may be performed 0- Decoherence may result from 
the interaction of the quantum system with a dissipa- 
tive environment 0] , dramatically affecting the dynamics 
of the quantum system, and its entanglement properties. 
Specifically, it was recently shown 0, Q that a (non inter- 
acting) two-qubit system can be completely disentangled 
in a finite time, a phenomenon dubbed "entanglement 
sudden death" . This was followed by a plethora of theo- 
retical studies of this phenomenon in various situations, 
most of them treating non-interacting qubits (i.e. a pair 
of qubits which interact with each other only in the me- 
diation of the environment), either with a Markovian or 
non-Markovian environment [1, 0, i, H, [13, [H El, El, [II ■ 
From the experimental side, direct measurements of en- 
tanglement have been performed [H, [l^. a nd the " en- 
tanglement sudden death" was observed pTf . 

The loss of entanglement seems to be a generic fea- 
ture of two-qubit systems [1|. Our goal is to study a 
simple system where entanglement sudden death may be 
avoided. For this aim we study a simple model of an 
interacting two-qubit system in the presence of a ther- 
mal dissipative bath. We model the qubits as interact- 
ing spin-i particles in a magnetic field [3] and use the 
Born-Markov approximation for the system-environment 
coupling. The environment is assumed to be ohmic and 
induce thermal transitions, and it only allows for sin- 

t spin-flip events (in similarity to spin-boson models 
). We calculate analytically the full dynamics of this 



system, with emphasis given to the different relaxation 
rates. We point out here that we use the term "relax- 
ation rates" loosely, to describe the time-scales of both 
processes which include energy changes (relaxation) and 
only coherence loss (decoherence). Both types of pro- 
cesses are inherently present in our calculation scheme. 
We show that although the coupling with the environ- 
ment is characterized by a single relaxation rate, differ- 
ent relaxation rates emerge for different coherent states. 
We study the effect of temperature and magnetic field on 
the different relaxation rates and find that they may be 
non-monotonous functions of temperature. 



As the main results of this paper, we demonstrate that 
as a function of the magnetic field the relaxation rate of 
the highly-entangled states abruptly changes from be- 
ing finite as the temperature vanishes to being exponen- 
tially small. This occurs when one of the states of the 
entangled pair is in a meta-stable state, and indicates 
that with a proper tuning of parameters the entangle- 
ment may survive very long times even in the presence 
of a dissipative environment. We demonstrate the long- 
lifeness of entanglement by calculating the concurrence 
of a specific entangled state, and show that applying a 
transverse magnetic field destroys this effect. 



Finally, we devise a simple way to account for non- 
Markovian effects when calculating relaxation rates. 
Studying the relaxation rates as a function of the environ- 
ment memory-time, we find that some relaxation rates 
exhibit a non-monotonic dependence on the memory- 
time, with a cusp-like resonance. For long memory-times, 
the different relaxation rates merge into a single rate. 
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II. METHOD 

Let us introduce our method for calculating the re- 
laxation times. We consider a quantum system, char- 
acterized by a time-independent Hamiltonian TL with 
iV energy levels k = 1,2,...N. The system dy- 
namics are given by the evolution of the density matrix 

Pit) = E 

kk' Pkk'{t)\k){k'\, where \k) are the eigenfunc- 
tions of the Hamiltonian. For the above choice of 
Hamiltonian, in the Markovian approximation the evolu- 
tion of the density matrix is given by a quantum master 
equation |20] {h= 1), 



(1) 



where Cp is a superoperator describing the dissipative 
dynamics. It is commonly takes the Lindblad form [2^ 



(2) 



(where {■, ■} are the anti-commutation relations), which 
ensures positivity of the density matrix The V- 

operators define the different relaxation processes in- 
duced by the environment. 

We now follow Ref. [25] and cast the den- 
sity matrix into a vector form, defining p — 
(pii, P22, ■■■,PNN, 3?pi2, 3pi2, ■■■,^piN,'^piN, ■■■)■ Here 
the first N elements account for occupation probabilities 
and the other elements describe coherence between the 
states in the statistical mixture. It is now a matter of re- 
arranging the master equation into a form p{t) = Mp(t), 
where now M is a matrix of dimension iV^ which in- 
cludes both the Hamiltonian and the dissipative part of 
the evolution. 

Due to the semi-group properties the Lindblad equa- 
tion, (at least) one of the eigenvalues of M is exactly zero 
[25| . The eigenvector corresponding to this eigenvalue is 
the steady-state of the system, i.e. the limit limj^oo p{t)- 
The other eigenvalues may have an imaginary part, but 
all of them have a negative real part, which correspond 
to the relaxation rate of the corresponding eigenvector. 
Thus, by calculating the eigenvalues of M one obtains 
the relaxation rates for all possible processes. For the 
most general initial condition, the smallest non- vanishing 
eigenvalue of M (which we call Ai) represents the longest 
relaxation rate. 



III. APPLICATION TO THE TWO-QUBIT 
SYSTEM 

A. Two-qubit system in a perpendicular field 

Next we consider the application of the above method 
to our model two-qubit system. The hamiltonian, with 
Ising interactions, can be written as 



1=1,2 



(3) 



where s^*) are the two qubit levels (i.e. spin=i par- 
ticles), J > describes the interactions and B is the 
external magnetic field. This choice of Hamiltonian is 
not only convenient (the Hamiltonian being very sim- 
ple) , but also represents several sugges tions for realistic, 
spin-based quantum computers |26l . |27| . 

For simplicity we start with magnetic field only in the 
z-direction, i.e. B = Bz with B > 0. Choosing as a 
basis the four states | tt), I Ti)J it),! ii) (which we 
number from 1 to 4, respectively), the Hamiltonian can 
be written (up to a constant energy shift) as 



n 



f -J -B 





V -J + S 



(4) 



In order to account for the environment, the spins are 
coupled to one that induces spin-flip processes, repre- 
sented by the F-operators in Eq. [51 Here we make two 
assumptions, namely (i) the spins are flipped one at a 
time (i.e. there is no direct relaxation from the | IT) 
state to the | [[) state, etc.) and (ii) the relaxation rate 
between two states is proportional to the Boltzman fac- 
tor of the corresponding energy difference between the 
states (these two assumptions on the form of the 
^-operators reflect the properties of the environment as 
described in the introduction). For example, the relax- 
ation operator from | ft) to | fj) (taking ks = 1) is 

El 



712 



7!2/'in)(TTi 

70 

cosh^^- 



( 2T ) 



J + B 
2T 



(5) 



where 70 is some typical relaxation rate which represents 
the strength of the qubit-environment coupling. For the 
reverse process the relaxation operator is 



^21 



721 



1/2 

721 



TT)(Ti 

70 



cosh ( 



J+B \ "^^P 
2T / 



J + B 
2T 



(6) 



so that the condition of detailed balance is maintained, 
i.e. 712/721 — exp (— ^^" ). Note that one can normal- 



ize the transition rate in different ways and still maintain 
detailed balance. Here we choose such a normalization 
that keeps all relaxation rates finite (even at T — 0) and 
preserves the B ^ —B symmetry (i.e., does not depend 
on the gauge of the Hamiltonian). However, all the qual- 
itative results presented in this paper equally apply for a 
different choice of the normalization of the F-operators. 

Once the form of the V^-operators is specified, it is 
now a matter of algebraic manipulation to obtain the 
M-matrix and its eigenvalues and eigenvectors. For the 
steady state we find 

p(^) = ^"'(lTT)(TTI + e"*^Ui)(iil 

-(J + B) 



+e" 



■(I n)(n I + UT)(iT I) , (7) 
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with Z being the partition function of the system, which 
results in a pure | IT) state (which is the ground state of 
the Hamiltonian) in the hmit of T ^ 0. 

For the above example, all the eigenvalues may be cal- 
culated analytically. For the lowest rate we find 



Ai -1 



sinh (y) 



cosh I 



cosh 



'2^ 



(8) 



which is 3-fold degenerate. Two states contain pi4 and 
/341 (i.e. a coherence between | IT) and | ||)), and the 
third is a mixture of all the diagonal elements pn, p22, P33 
and P44. Inspection of Ai in the limit T — > shows that 



lim Ai 



0, B <.J 

lol"^. B = J 
70, B> J 



(9) 



This means that for B < J the relaxation time from the 
coherent pi4 state diverges, i.e. the system never reaches 
its thermal, disentangled, ground state. However, for 
S > J it relaxes to the ground state in a time scale 
T ~ 7(7^- This can be easily explained from the following 
considerations. Note that the | H) state can only relax 
into one of the degenerate middle states, | ti) or | 11). 
For < B < J, its energy is negative (but higher than the 
ground state energy), and therefore the relaxation rate is 
exponentially small. Put it differently, the | H) state is 
a meta-stable state which requires an exponentially rare 
correlated event to escape from. For B > J, however, 
the energy of | H) is no longer negative, and hence it 
can relax into the ground state by a cascade relaxation 
through the middle states. 

While all the other eigenvalues are available analyt- 
ically, writing them in full form is cumbersome, and 
thus we present them graphically. In Fig. [T] we plot 
the eigenvalues (corresponding to the inverse relaxation 
rates for the different states), as a function of temper- 
ature for magnetic field, B/J = 0.9. In Fig. [TJa) we 
plot a wide temperature range, and we zoom in on the 
low-temperature range in Fig. [Ijb). For each relaxation 
rate the corresponding elements of the density matrix are 
marked. Fig.[l]shows two interesting features: (i) the re- 
laxation rates are nonmonotonic in temperature, and (ii) 
they break into groups, with only two possible time scales 
( 7(7^ and 270"^) at T -> 0. 

For comparison, in Fig. [2] we plot the same for a mag- 
netic field B/J = 1.1. We now find that the relaxation 
rates are broken into four groups, i.e. additional time- 
scales appear. Interestingly, only the /923-state preserves 
its T — > limit as the B = J point is crossed, while all 
other time-scales exhibit a discontinuous change. 



B. Concurrence 

In order to demonstrate the above effect on the entan- 
glement of the two-qubit system, we calculate the con- 
currence, which is a direct measure of the entanglement 
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FIG. 1: (color online) different relaxation rates (at B/J = 
0.9) as a function of temperature for (a) a large tempera- 
ture scale, and (b) zoomed in on the low-temperature regime. 
For each rate, the corresponding density matrix element is 
pointed. 

^2§\ . For the two-qubit system (in the basis chosen above) 
it is defined as C = max(0, — — — y^), 
where e,; are the eigenvalues of the matrix C, defined by 



J. {'2) J. z 

C — p{s)j' ® Sy )p*{sy'^' (g) s\f'). Clearly, the dynamics of 
the concurrence depend on the initial condition. For this 
example we choose the initial density matrix 



„(1) ^ J2). 



p(0) = 



/ pii pi4 

P22 

P33 

V P41 P44 



(10) 



for which one can easily calculate the concurrence, 
C{t) = max(0, 2(^pi4p4i - ^/piipii). The diagonal el- 
ements are given by Eq. ([7]) and pi4 — pn = 1/2. By 
making this choice we start with a highly entangled state 
(C ~ 1), but the dynamics are very easy to calculate as 
the diagonal elements do not change at all, and the off- 
diagonal ones decay with the rate given by Eq. ([9]). We 
point out that one can start with any initial diagonal 
elements and obtain results similar to those presented 
below, since the diagonal elements quickly relax to the 
steady-state (Eq. (O) and do not contribute to the con- 
currence time-evolution. 

In Fig. [3] we plot the concurrence as a function of tem- 
perature and time for magnetic field values B/J = 0.9 
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FIG. 2: (color online) Same as in Fig. [T]but for B/J = 1.1 



(upper panel) and B/J ~ 1.1 (lower panel). The mesh 
corresponds to finite C, while in the blank regions C — 0. 
For B > J we find that the " entanglement sudden death" 
is present at all temperatures. However, for B < J it be- 
comes exponentially suppressed at lower temperatures: 
the concurrence practically remains finite for all times. 

In order to understand this behavior, we note that the 
formula for the concurrence describes a competition be- 
tween the off-diagonal elements of the density matrix, 
which contribute to the entanglement, and the diago- 
nal elements, which "disentangle" the state. At strictly 
T = 0, the diagonal elements vanish but the off-diagonals 
survive indefinitely, giving rise to an entangled state. For 
very low temperature, while in the strict t — s- 00 limit the 
system becomes disentangled, this time is exponentially 
long. 



C. Transverse field 

Let us consider the effect of an additional transverse 
field B — i?2.x (the results are identical to an additional 
field in the y-direction). A transverse field can either 
represent an inherent interaction between the two qubit 
states, or an actual field (in the qubit relevant Hilbert 
space), which is used to perform quantum operations. 
Since both these ingredients appear in any implementa- 
tion of a physical qubit, it is important to study their 
effect on the relaxation time-scales. 




FIG. 3: Concurrence C as a function of temperature and 
time for magnetic field values B/J = 0.9 (upper panel) and 
B/J = 1.1 (lower panel). The Mesh corresponds to finite 
C, while in the blank regions C — 0. For B < J at low 
temperatures the concurrence practically never decays, while 
for B > J sudden death of entanglement appears for every 
temperature. 



We thus repeat the above procedure of constructing 
and diagonalizing the M-matrix with an additional term. 
In Fig. m we plot the inverse relaxation rates at T = 
0,B = 0.9 as a function of B^- It is found that the infinite 
relaxation time becomes finite (i.e. there is no more a 
meta-stable state) and that the degeneracy is lifted. This 
means that applying a transverse field might give rise to 
entanglement " sudden-death" . Interestingly, one can see 
that the time-scales are not monotonic functions of B^- 



The conclusion arising from the above calculation is 
two-fold, (i) in a realistic qubit, coupling between the 
qubit states should be maximally inhibited, to allow for 
longer entanglement life-time, and (ii) one has to take 
into account the fact that performing quantum opera- 
tions on the qubits will result in faster decay of entangle- 
ment. 
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FIG. 4: (color online) Relaxation rates (at T = 0, B / J = 0.9) 
as a function of perpendicular magnetic field Bx- 

IV. NON-MARKOVIAN EFFECTS 

In the Markovian approximation, the evolution of the 
system does not depend on its history. This approxi- 
mation applies when the correlation times of the thermal 
bath are much smaller than any timescale associated with 
the system under consideration. However, there are cases 
where the bath reacts to the dynamics of the system over 
a certain memory time, tm, and the Markov approxima- 
tion is thus not valid. In such instances, the evolution 
of the system would depend on its history [23|- Non- 
Markovian effects on the entanglement dynamics have 
been extensively studied in recent years [ll|, [H, [l^ . 

In order to include non-Markovian effects, one has to 
include the history of the system. In the simplest ap- 
proximation (23 | , this adds up to a form of the quantum 
master equation 

p(t) = -i[H,p]+ / K{t-t')Cp{t')dt' , (11) 
Jo 

where the memory kernel K{t) defines the response of the 
bath to the history of the system. Note that not every 
form of the memory kernel is possible, as positivity of the 
density matrix may be lost [29| . 

One can now repeat the procedure described above, 
and rewrite Eq. (jlip in a vector form, which yields the 
integro-differential vector equation 

'p{t) = / M{t ~ t')p{t')dt' . (12) 
Jo 

The time-dependence of the M-matrix is such that for 
elements derived from the Lindblad super-operator one 
attaches the kernel K(t—t'), and for the elements derived 
from the Hamiltonian one attaches a (5-function, 5{t — t'). 
One can now Laplace-transform Eq. ^\2\ and obtain 
an algebraic vector-equation sp{s) — p{Q) ~ M{s)p{s), 
where s is the (complex) Laplace variable. Note that 
the Laplace transform M(s) has a simple form, as the 



Hamiltonian elements are multiplied by unity (which is 
the transform for the J-function) and the Lindbladian 
elements are multiplied by the Laplace transform oiK{t). 
The formal solution of the above equation is thus given 

by 

/ioc 
e"*(sl- M(s))-V(0)ds , (13) 
-ioQ 

where I is the unit matrix. From Eq. (jl3p it can be seen 
that if the secular equation dot (si — M{s)) has solutions, 
the real part of these solutions defines a relaxation time- 
scale (via a Cauchy-like integration over poles) . One can 
thus obtain the relaxation times from a numerical solu- 
tion of the secular equation, without the need to solve 
the full non-Markovian dynamics. 

In the simplest approximation [2^, [s^, HH, [12] the mem- 
ory kernel is given by an exponential form, 

if(t)=T-/exp(-|i|/rM) , (14) 

where tm is the memory-time. The Laplace transform 
is thus K{s) = ^ , which yields a polynomial secular 
equation. Note tliat in the limit tm the Markovian 
limit is exactly obtained. We have solved this equation 
numerically, and found that it always has N — 1 solu- 
tions with negative real part and a single solution with 
3 = 0, corresponding to the steady-state. This means 
that for an exponential memory-kernel, one can identify 
different processes which have different relaxation rates. 
We note that if one takes a more complicated memory 
kernel (say a power-law), then the secular equation be- 
comes transcendental, with no simple poles. In that kind 
of environment, one cannot simply attach different time- 
scales to different processes, and the full dynamics of the 
system must be calculated [33j . 

The above scheme can now be applied to our two-qubit 
model. In Fig. [5] the inverse time-scales are plotted as a 
function of the memory-time tm for B/J = 0.9, 1.1. We 
find that the time-scales have a resonant-structure as a 
function of tm, and in fact exhibit a cusp at the reso- 
nance. In addition, at relatively large tm, all the differ- 
ent time-scales merge into a single group, i.e. there is no 
longer a separation of the time-scale for the relaxation 
processes of different coherent states. 

The long memory-time behavior may be understood 
by the fact that in this case the bath memory dominates 
the relaxation processes, giving rise to a single time-scale. 
The longer the memory-time, the more weight is given 
in the relaxation process to states far from the steady- 
state, and hence the relaxation rate diminishes. However, 
the rise in relaxation rate at small memory-times is a 
surprising effect, which comes about due to the complex 
nature of the interaction between the two-qubit system 
and the environment. 
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FIG. 5: (color online) Relaxation rates (at T — 0,B/J — 
0.9, 0.1) as a function of the memory-time tm- Note the reso- 
nant structure of the time-scales, and the bunching of different 
time-scales at large r. 

V. SUMMARY 

In this paper we have studied the relaxation dynamics 
of a model interacting two-qubit system in the presence 
of an ohmic thermal environment. The qubits were mod- 
eled as spin-i particles with spin-spin coupling and in 
the presence of a magnetic field. The environment was 
limited to induce only single spin- flip events. Within this 
model we analytically calculated the relaxation rates of 
different processes. Our main result is that disentangle- 
ment may be critically slowed down in the T — > limit 
by varying the magnetic field, and entanglement sudden 
death may be completely avoided (or at least exponen- 



tially suppressed for low temperatures). This was ex- 
plicitly shown by calculating the concurrence dynamics 
of a highly entangled state for different magnetic fields. 
We have also shown that a transverse magnetic field may 
destroy this effect. 

In addition, we have introduced a simple way to include 
non-Markovian effects in the calculation of the relaxation 
rates. We have shown that the different relaxation rates 
exhibit an interesting non-monotonic dependence on the 
environment memory-time, with a cusp-like resonance. 
For long memory-times, we have found that the different 
time-scales merge into a single time-scale. 

In order to experimentally verify the results presented 
here, one needs an experiment with a well-controlled two- 
qubit system, where the energy-levels can be controlled. 
A promising candidate is a system composed of two cou- 
pled two-quantum-dot qubits (33 | , where the qubit levels 
may be controlled by external gates, and a high level of 
control of the qubit states has been already demonstrated 
experimentally [ssj . 
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